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Abstract. Graph kernels methods are based on an implicit embedding 
of graphs within a vector space of large dimension. This implicit embed- 
ding allows to apply to graphs methods which where until recently solely 
reserved to numerical data. Within the shape classification framework, 
graphs are often produced by a skeletonization step which is sensitive 
to noise. We propose in this paper to integrate the robustness to struc- 
tural noise by using a kernel based on a bag of path where each path 
is associated to a hierarchy encoding successive simplifications of the 
path. Several experiments prove the robustness and the flexibility of our 
approach compared to alternative shape classification methods. 
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1 Introduction 

The skeleton of a 2D shape is defined as the location of the singularities of the 
signed distance function to the border of the shape. This structure has several 
interesting properties: it is thin, homotopic to the shape, invariant under rigid 
transformations of the plane and most importantly it has a natural interpretation 
as a graph. The representation of a shape by a skeletal (or shock) graph has 
become popular owing the good properties of this representation inherited from 
the properties of the skeleton. However, beside all this good properties, the 
skeletonization is not continuous and small perturbations of the boundary insert 
structural noise within the graph encoding the shape. 

Several graph based methods have been proposed to compute a distance be- 
tween shapes robust to such a structural noise. Sharvit et al. [1] propose a graph 
matching method based on a graduated assignment algorithm. Siddiqi [2] pro- 
poses to transform the shock graph into a tree and then applies a tree matching 
algorithm. Pcllilo [3] uses the same tree representation but transforms the tree 
matching problem into a maximal clique problem within a specific association 
graph. 



This work is performed in close collaboration with the laboratory Cyceron and is 
supported by the CNRS and the region Basse-Normandie. 
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All the above graph methods operate directly on the space of graphs which 
contains almost no mathematical structure. This lack of mathematical structure 
forbids the use of basic statistical tools such as the mean or the variance. Graph 
kernels provide an elegant solution to this problem. Using appropriate kernels, 
graphs can be mapped either explicitly or implicitly into a vector space whose 
dot product corresponds to the kernel function. All the "natural" operations 
on a set of graphs which were not defined in the original graph space are now 
possible into this transformed vector space. In particular, graph kernels may be 
combined with the kcrncliscd version of robust classification algorithms such as 
the Support Vector Machine (SVM). 

A Graph kernel used within the shape representation framework should 
take into account the structural noise induced by the skeletonization process. 
Bunke [4] proposes to combine edit distance and graph kernels by using a set 
of n prototype graphs {c/i, . . . , g n }. Given a graph edit distance .), Bunke 
associates to each graph g the vector 4>(g) = (d(g,gi) . . . ,d(g,g n )). The kernel 
fc(<?l,<72) between the two graphs g\ and 52 is then defined as the dot product 
< 0(si),0(sa) >. 

Neuhaus [5] proposes a similar idea by defining for a prototype graph go, the 
kernel: k go (g,g') = \ (d 2 (g,g ) + d 2 {g ,g') - d 2 {g,g')), where d(., .) denotes the 
graph edit distance. Several graph prototypes may be incorporated by summing 
or multiplying such kernels. Using both Neuhaus [5] and Bunke [4] kernels two 
close graphs should have close edit distance to the different graph prototypes. 
The metric induced by such graph kernels is thus relative both to the weights 
used to define the edit distance and to the graph prototypes. This explicit use 
of prototype graphs may appear as artificial in some application. Moreover, the 
definite positive property of these kernels may not in general be guaranteed. 

Suard [6] proposes to use the notion of bag of paths of finite length for shape 
matching. This method associates to each graph all its paths whose length is 
lower than a given threshold. The basic idea of this approach is that two close 
shapes should share a large amount of paths. A kernel between these sets should 
thus reflect this proximity. However, small perturbations may drastically reduce 
the number of common paths between two shapes (Section 3.2). Moreover, the 
straightforward definition of a kernel between set of paths does not lead to a 
definite positive kernel (Section 2.1). 

This paper proposes a new definite positive kernel between set of paths which 
takes into account the structural noise induced by the skeletonization process. 
We first present in Section 2 the bag of paths approach for shape similarity. Our 
contributions to this field are then presented in Section 3. The effectiveness of 
our method is demonstrated through experiments in Section 4. 

2 Kernels on bag of paths 

Let us consider a graph G = (V, E) where V denotes the set of vertices and 
E C V x V the set of edges. As mentioned in Section 1 a bag of paths P of 
length s associated to G contains all the paths of G of length lower than s. We 
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denote by \P\ the number of paths inside P. Let us denote by K pat h a generic 
kernel between paths. Given two graphs G\ and Gi and two paths h\ G P\ 
and hi € P2 of respectively G\ and G2, K pat h(hi, hi) may be interpreted as a 
measure of similarity between h± and /12 and thus as a local measure of similarity 
between these two graphs. The aim of a kernel between bags of paths consits to 
agregate all these local measures between pairs of paths into a global similarity 
measure between the two graphs. 

2.1 The max kernel 

This first method, proposed by Suard [6], uses the kernel K pat h as a measure of 
similarity and computes for each path hi 6 Pi the similarity with its closest path 
in Pi (max^ e p 2 K pat hih\, hj)). A first global measure of similarity between Pi 
and Pi is then defined as: 



The function K max {G\^Gi) is however not symmetric according to G\ and 
Gi. Suard obtains a symmetric function interpreted as a graph kernel by taking 
the mean of K max (Gi,Gi) and K max {G 2 ,Gi): 



This kernel is not positive definite in general. However as shown by Haasdonk 
[7] , SVM with indefinite kernels have in some cases a geometrical interpretation 
as the maximization of distances between convex hulls. Moreover, experiments 
(section 4, and [6]) show that this kernel usually leads to valuable results. 

2.2 The matching kernel 

The non definite positiveness of the kernel K max is mainly due to the max oper- 
ator. Suard [6] proposes to replace the kernel K pat h by a kernel which decreases 
abruptly when the two paths are different. The resulting kernel is defined as: 

Kmatching(Gi, G 2 ) = K ma tching(Pl7 -f*2 ) = (3) 



where d pa th is the distance associated to the kernel K pat h and defined by: 
d%ath(hi,h 2 ) = K pa th{hi,hi) + K pa th(h 2 , h 2 ) - 2K pat h{hi,h 2 ). 

The resulting function defines a definite positive kernel. This kernel relies 
on the assumption that using a small value of er, the couple of paths with the 
smallest distance will predominate the others in equation 3. This kernel may 
thus lead to erroneous results if the distance are of the same order of magnitude 
than a or if several couples of paths have nearly similar distances. 




(1) 



K max (Gi,Gi) — \ K m ax{Gi,G 2 ) + K max (Gi,Gi) 



(2) 
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(a) Sets on the unit (b) original (c) edge protru- (d) node 
sphere sion insertion 



Fig. 1. Separating two sets using one-class SVM (a). The symbols (wi, pi) and (w2, P2) 
denote the parameters of the two hyperplanes which are represented by dashed lines. 
Influence of small perturbations on the bag of paths ((b), (c) and (d)) 



2.3 The change detection kernel 

Desobry [8] proposed a general approach for the comparison of two sets which 
has straightforward applications in the design of a kernel between bags (sets) 
of paths. Desobry models the two sets as the observation of two sets of random 
variables in a feature space and proposes to estimate a distance between the two 
distributions without explicitly building the pdf of the two sets. 

The feature space considered by Desobry is based on the normalised kernel 
(K(h,h') = K pat h(h, h')/y/ (K pa th(h, h)K pat h(h', h'))). Using such a kernel we 
have = K(h, h) = 1 for any path. The image in the feature space of our set 

of paths lies thus on an hypersphere of radius 1 centered at the origin (Fig. 1). 
Desobry defines a region on this sphere by using a single class v-SVM. This 
region corresponds to the density support estimate of the unknown pdf of the 
set of paths [8] . 

Using Desobry's method, two set of vectors are thus map onto two regions 
of the unit sphere and the distance between the two regions corresponds to a 
distance between the two sets. Several kernels based on this mapping have been 
proposed: 

1. Desobry proposed [8] to define the distance between the two spherical arcs as 

/ »iK 12 «j \ 

aiXCOS I j. rrj. n- I 

a contrast measure defined by: d 2 Desob (P 1 ,P 2 ) = M^jngay _ 

This distance is connected to the Fisher ratio (see [8, Section IV]). However, 
the definite positiveness of the Gaussian RBF kernel based on this distance 
remains to be shown. 

2. Suard [6] proposed the following kernel: Ksi ia rd(Gi,G 2 ) = K Su ard(Pi, P2) = 
PiP^J2f H ePi 52h j eP 2 ai - i K path(hi,h.j)a2j with wi = (ai,i, . . . , 0-1,^1) and 
w 2 = (a2,i, ■ ■ .,a 2 ,|f(j|). 

This kernel is definite positive, but does not correspond to any straightfor- 
ward geometric interpretation. 
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2.4 Path kernel 

All the kernels between bags of paths denned in Section 2 are based on a generic 
kernel K pat h between paths. A kernel between two paths h\ = (v±, . . . , v n ) and 
b! = (v[, ... , v' p ) is classically [9] built by considering each path as a sequence of 
nodes and a sequence of edges. This kernel denoted K c i assic is then defined as 
if both paths have not the same size and as follows otherwise: 

\h\ 

K c ia SS ic(h, ti) = K v (ip(vi), ip(v' 1 ))'jJ l _K e (ip(e Vi _ 1Vi ), ip(e v <. _ iV 0)K v (ip(vi),ip(vi)) 

i=1 

(4) 

where ip{v) and 2p(e) denote respectively the vectors of features associated to 
the node v and the edge e. The terms K v and K e denote two kernels between 
respectively nodes and edge's features. For the sake of simplicity, we have used 
Gaussian RBF kernels between the attributes of nodes and edges (Section 4). 

3 Hierarchical kernels 

Since the main focus of this paper is a new kernel method for shape classification, 
the construction of skeletal graphs from shapes has been adressed using classical 
methods. We first build a skeleton using the method proposed by Siddiqi [2]. 
However the graph we build from the skeleton does not correspond to the shock 
graph proposed by Siddiqi. Indeed, this graph provides a precise description 
of the shape but remains sensitive to small perturbations of the boundary. We 
rather use the construction scheme proposed by Suard [6] and Ruberto [10] which 
consists to select as node all the pixels of the skeleton which correspond to end 
points or junctions. These nodes are then connected by edges, each edge being 
associated to one branch of the skeleton. Given a skeletal graph G we valuate 
each of its edge by an additive weight measure and we consider the maximal 
spanning tree T of G. The bag of path associated to G is built on the tree T. 
Note that, the skeletonization being homotopic we have G = T if the 2D shape 
does not contain any hole. 

3.1 Bag of path kernel 

None of the bag of path kernels proposed by Desobry or Suard (Section 2) is 
both definite positive and provides a clear geometrical interpretation. We thus 
propose a new kernel based on the following distance: 

d change (Pi,P2) = arccos ( ^^^ J • (5) 

This distance corresponds to the angle a between the two mean vectors w\ and 
W2 of each region (Fig. 1). Such an angle may be interpreted as the geodesic 
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distance between two points on the sphere and has thus a clear geometrical 
interpretation. Based on this distance we use the Gaussian RBF kernel: 

K change {Gl, G 2 ) = K change (P 1 ,P 2 ) = CXp ^ change^ 1> gA ^ ^ 

This kernel is definite positive since the normalized scalar product is positive 
definite and arccos is bijective on [0, 1]. The Gaussian RBF kernel based on this 
distance is thus definite positive (see [11] for further details). 



3.2 Hierarchical kernel between paths 

A mentioned in Section 1, the use of kernels between bags of paths within the 
shape matching framework relies on the assumption that the graphs associated 
to two similar shapes share a large amount of similar paths. This assumption 
is partially false since a small amount of structural noise may have important 
consequences on the set of paths. Let us for example, consider the small defor- 
mation of the square (Fig. 1(b)) represented on Fig. 1(c). This small deformation 
transforms the central node in Fig. 1(b) into an edge (Fig. 1(c)). Consequently 
graphs associated to these two shapes only share two paths of length 2 (the ones 
which connect the two corners on the left and right sides). In the same way, 
a small perturbation of the boundary of the shape may add branches to the 
skeleton(Fig. 1(d)). Such additional branches i) split existing edges into two sub 
edges by adding a node and ii) increase the size of the bag of path either by 
adding new paths or by adding edges within existing paths. 

The influence of small perturbations of the shape onto an existing set of paths 
may thus be modeled by node and edge insertions along these paths. In order 
to get a path kernel robust against structural noise we associate to each path a 
sequence of successively reduced paths, thus forming a hierarchy of paths. Our 
implicit assumption is that, if a path has been elongated by structural noise one 
of its reduced version should corresponds to the original path. 

The reduction of a path is performed either by node removal or edge con- 
traction along the path. Such a set of reduction operations is compatible with 
the taxinomy of topological transition of the skeleton compiled by Giblin and 
Kimia [12]. Note that, since all vertices have a degree lower than 2 along the path 
these operations are well defined. In order to select both the type of operation 
and the node or the edge to respectively remove or contract we have to associate 
a weight to each node and edge which reflects its importance according to the 
considered path and the whole graph. 

Let us consider a skeletal graph G, its associated maximal spanning tree T 
and a path h = (vi, . . . , v n ) within T. We valuate each operation on h as follows: 

Node removal : Let us denote by Vi, i <E {2, . . . , n — 1} the removed node of 
the path h. The node u, has a degree greater than 2 in T by construction. 
Our basic idea consists to valuate the importance of Vi by the total weight 
of the additional branches which justify its existence within the path h. For 
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(a) original skeleton (b) reduction of node 2 (c) reduction of edge e 2 .;j 



Fig. 2. Different reductions of a path (in gray) within a skeletal tree. 



each neighbor v of Wj not equal to nor Vi+i we compute the weight W(v) 
defined as the addition of the weight of the tree rooted on v in T — {e ViV } 
and the weight of e ViV . This tree is unique since T is a tree. The weight of 
the node Vi (and the cost of its removal) is then defined as the sum of weight 
W(v) for all neighbors v of Vi (excluding Vi—i and vi). 

After the removal of this node the edges e Vi _ lVi and e ViVi+1 are concatenated 
into a single edge in the new path h! . The weight of this new edge is defined 
as the sum of the weight of the edges e Vi _ lVi , e ViVi+1 and the weight of the 
node Vi (Fig. 2(a) and (b)). 
Edge contraction : The cost of an edge contraction is measured by the rele- 
vance of the edge which is encoded by its weight. Let us denote by e ViVi+1 , 
i < n the contracted edge of the path h = (v\, . . . ,v n ). In order to preserve 
the total weight of the tree after the contraction, the weight of the edge 
£viv i+1 is equally distributed among the edges of T incident to Vi and v i+1 : 

Ve e t(vi) U t(v i+1 ) - {e ViVi+1 } w'{e) = w(e) + 



d(vi) + d(v i+ i) - 2 

where i(v) and d(v) denote respectively the set of edges incident to v and 
the cardinal of this set (the vertex's degree). The symbol w(e) denotes the 
weight of the edge e. 

For example, the contraction of the edge e2,3 in Fig. 2(a) corresponds to a 
cost of w(e2,3) = -5. The contraction of this edge induces the incrementation 
of the edge's weights w(e2,i),w(e2,6),w{s3,4:) by -5/3 « .16 

Any additive measure encoding the relevance of a branch of the skeleton may 
be used as a weight. We choose to use the measure defined by Torsello [13] which 
associates to each branch of the skeleton (an thus to each edge) the length of 
the boundaries which contributed to the creation of this branch. Such a measure 
initially defined for each pixel of the skeleton is trivially additive. 

Let us denote by k the function which applies the cheapest operation on a 
path. The successive applications of the function k associate to each path h a 
sequence of reduced paths [h, n(h), . . . , n D {h)) where D denotes the maximal 
number of reductions. Using K c i ass i C for the path comparison, we define the 
kernel K ec tit as the mean value of kernels between reduced paths of equal length. 
Given two paths h and h', this kernel is thus equal to if \\h\ — \h'\\ > D. Indeed, 
in this case the maximal reduction of the longucst path remains longuer than 
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the shortest one. Otherwise, \\h\ — \h'\\ < D, and K e dit{h, h') is defined as: 



This kernel is proportional (by a factor D + 1) to a sum of i?-convolution 
kernels [14, Lemma 1] and is thus definite positive. 

Since K c i ass i c is equal to for paths of different lengths, K e dit is indeed equal 
to a sum of kernels between reduced paths of equal length. For example, given 
two paths h and hi whose respective length is equal to 4 and 3 we have for D = 2: 



4 Experiments 

We used the following features for our experiments: Each node is weighted by its 
distance to the gravity center of the shape and each edge is assocated to a vector 
of two features: The first feature corresponds to the edge's weight (section 3.2). 
The second feature is the angle between the straight line passing through the 
two nodes of the edge and the principal axis of the shape. These experiments 
are based on the LEMS [15] database which consists of 99 objects divided into 
9 classes. 

We defined three kernels for these experiments: The kernel K m ax. classic based 
on a conjoint use of the kernels K max (equation 2) and K c i ass i C (equation 4) has 
been introduced by Suard [6]. The kernel K c i iange . c iassic based on a conjoint use of 
the kernels K c hange (equation 6) and K c i ass i C allows to evaluate the performances 
of the kernel K c hange compared to the kernel K max . Finally, the kernel K new is 
based on a conjoint used of the two kernels K e dit and Kchange proposed in this 
paper. The kernel Kciassic is defined by the two parameters <j e dge and (J vertex 
respectively used by the Gaussian RBF kernels on edges and vertices. The kernel 
Kmax,classic does not require additional parameters while K change , c iassic is based 
on a f-SVM and requires thus the parameter v. It additionally requires the 
parameter uj:^™ used by the RBF kernel in equation 6. The kernel K e( nt 
requires the two parameters <j e d ge and a ver t e x used by K c i ass i C together with the 
maximal number of edition (D). Finally, the kernel K new requires as K c hange 
the two additional parameters v and (J change (equation 6). These parameters have 
been fixed to the following values in the experiments described below: D = 2, 



(Tedge = Vertex = 0.1, v = 0.9, a^ nge = 0-3 and af^ g % = 1.0. The parameters 



a e dge and a ver tex are common to all kernels. The remaining parameters have 
been been set in order to maximize the performances of each kernel on the 
experiments below. 

Our first experiment compares the distance induced by each kernel k and 
defined as d 2 (x,x') = k{x, x)+k(x l ,x')—2k(x,x'). The mean number of matches 
for each class is defined as follows: For each shape of the selected class we sort 
all the shapes of the database according to their distances to the selected shape 




fc=0 1=0 





OS SIC 



(K 2 (h),K(h>))] 
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^> - s — v (2) if 
ft ft ft ft ft (3) 4, 

(a) hands, tools, dudes 



^ t 1 = 7 

(b) sorted distances to the hand 



Fig. 3. Five representative shapes of the classes hands, tools and dudes of the LEMS 
database (a), and (b) the 10 closest shapes from an hand using the distances induced 
by the kernels K maXyC i assic (1), K ch 

ange. classic 

(2) and Knew (3). 



using an ascending order. The number of good matches of the input shape is then 
defined as the number of shapes ranked before the first shape which belongs to 
a different class than the selected one. For example, the 10 nearest neighbors of 
a hand sorted in an ascending order are represented in Fig. 3(b), the number 
of good matches of each shape is indicated on the right of the figure. Note 
that the greater number of good match being obtained for the kernel K new . 
The mean number of good matches of a class is defined as the mean value of 
the number of good matches for each shape of the class. The different values 
represented in Tab. 1(a) represent the mean values of these number of good 
matches for the classes: hands, tools and dudes (Fig. 3(a)). As indicated by 
Tab. 1(a), the kernel K maX:C i assic provides stable results but is sensitive to slight 
perturbations of the shapes as the ones of the class dudes and cannot handle the 
severe modifications of the hands. The kernel K c h a nge,ciassic leads to roughly 
similar results on the different classes. Though not presented here, the kernel 
K matching, classic (equation 3) gives worst results than the others kernels. This 
result may be explained by the drawbacks of this kernel (Section 2.2). The kernel 
K new always provides the best results with a good robustness to perturbation 
on dudes and hands. 

Our second experiment evaluates performances of each kernel within a clas- 
sification framework. To this end, we have trained a SVM on 5 shapes of each of 
the three classes: dudes, hands, tools on one side and one model of each of the 
6 remaining classes on the other side. The SVM margin parameter was selected 
in order to maximize the number of true positive while having no false positive. 
Tab. 1(b) shows the number of well classified shapes for each class. The kernel 





Hands 


Tools 


Dudes 


Hands 


Tools 


Dudes 


K max , classic 


4.81 


6.18 


6.36 K max , classic 


7 


11 


10 


K change, classic 


5.27 


5.45 


6.36 K change, classic 


7 


10 


10 


Knew 


7.09 


9.82 


6.36 Knew 


9 


11 


11 



(a) Mean number of good (b) Number of recognized shapes in 

matches. one class. 

Table 1. Kernels evaluation based on distance (a) and classification (b) criteria. 
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Knew gives the best performances especially for the hands where the two missing 
shapes are the more perturbed ones. The two others kernels present good results 
and arc competitive when shapes arc not strongly deformed. This experiment 
confirms the robustness of our kernel against perturbed shapes. 

5 Conclusion 

The bag of path approach is based on a decomposition of the complex graph 
structure into a set of linear objects (paths). Such an approach benefits of recent 
advances in both string and vectors kernels. Our graph kernel based on a hier- 
archy of paths is more stable to small perturbations of the shapes than kernels 
based solely on a bag of paths. Our notion of path's hierarchy is related to the 
graph edit distance through the successive rewritings of a path. Our kernel is 
thus related to the ones introduced by Ncuhaus and Bunkc. 
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Preface 



This textbook is intended for use by students of physics, physical chemistry, 
and theoretical chemistry. The reader is presumed to have a basic knowledge 
of atomic and quantum physics at the level provided, for example, by the first 
few chapters in our book The Physics of Atoms and Quanta. The student of 
physics will find here material which should be included in the basic education 
of every physicist. This book should furthermore allow students to acquire an 
appreciation of the breadth and variety within the field of molecular physics and 
its future as a fascinating area of research. 

For the student of chemistry, the concepts introduced in this book will provide 
a theoretical framework for that entire field of study. With the help of these con- 
cepts, it is at least in principle possible to reduce the enormous body of empirical 
chemical knowledge to a few basic principles: those of quantum mechanics. In 
addition, modern physical methods whose fundamentals are introduced here are 
becoming increasingly important in chemistry and now represent indispensable 
tools for the chemist. As examples, we might mention the structural analysis of 
complex organic compounds, spectroscopic investigation of very rapid reaction 
processes or, as a practical application, the remote detection of pollutants in the 
air. 
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Abstract. The abstract should summarize the contents of the paper 
using at least 70 and at most 150 words. It will be set in 9-point font 
size and be inset 1.0 cm from the right and left margins. There will be 
two blank lines before and after the Abstract. . . . 

1 Fixed-Period Problems: The Sublinear Case 

With this chapter, the preliminaries are over, and we begin the search for periodic 
solutions to Hamiltonian systems. All this will be done in the convex case; that 
is, we shall study the boundary-value problem 

x = JH'(t,x) 
x(0) = x(T) 

with H(t, •) a convex function of x, going to +00 when ||x|| — ► 00. 
1.1 Autonomous Systems 

In this section, we will consider the case when the Hamiltonian H(x) is au- 
tonomous. For the sake of simplicity, we shall also assume that it is C x . 

We shall first consider the question of nontriviality, within the general frame- 
work of (^4oo, i?oo)-subquadratic Hamiltonians. In the second subsection, we shall 
look into the special case when H is (0, 6 oc )-subquadratic, and we shall try to 
derive additional information. 

The General Case: Nontriviality. We assume that H is {A^, _B oc )-sub- 
quadratic at infinity, for some constant symmetric matrices A^ and B^, with 
Boo — Aoo positive definite. Set: 



7 : = smallest eigenvalue of — A, 



00 



(1) 



A : = largest negative eigenvalue of J— + A, 



dt 



OO 



(2) 



Theorem 1 tells us that if A + 7 < 0, the boundary- value problem: 

x = JH'(x) 
x(0) = x{T) 



(3) 



has at least one solution x, which is found by minimizing the dual action func- 
tional: 

" "1 



(^jVu) +i\T(-u) 



dt 



(4) 



on the range of A, which is a subspace R(A) L with finite codimension. Here 



N(x) :=H(x)--(A oc x,x) 

is a convex function, and 

N(x) < i ((Boo -Aoo)x,x) +c \fx 
Proposition 1. Assume H'(0) = and H(0) = 0. Set: 

6 := liminf 2N{x) ||x|f 2 . 

7/7 < —A < 5, the solution u is non-zero: 

x{t) ^0 Vi . 



(5) 



(6) 



(7) 



(8) 



Proof. Condition (7) means that, for every 8' > 5, there is some e > such that 

(9) 



Ml ^^(z)^ ||x|| 2 



It is an exercise in convex analysis, into which we shall not go, to show that 
this implies that there is an i] > such that 



f\\x\\<V^N*(y)<^j\\y\\ 2 



(10) 



Fig. 1. This is the caption of the figure displaying a white eagle and a white horse on 
a snow field 
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Since u\ is a smooth function, we will have H/wiiH^ < i] for h small enough, 
and inequality (10) will hold, yielding thereby: 



ip(hut) < \\ \\ Ul \\l + ^-^ ||ui|| 2 



(11) 



If we choose 5' close enough to S, the quantity ( j + jr) will be negative, and 
we end up with 

ij){hu{) < for h ^ small . (12) 

On the other hand, we check directly that tp(0) = 0. This shows that cannot 
be a minimizer of tjj, not even a local one. So u ^ and u ^ A~ 1 (0) = 0. □ 

Corollary 1. Assume H is C 2 and (doo, &oo) -subquadratic at infinity. Let £i, 
...,£n be the equilibria, that is, the solutions of H'(£) = 0. Denote by uik the 
smallest eigenvalue of H" (£fc), and set: 



uo := Min {u)\, . . . , w*;} 



T 



(13) 
(14) 



then minimization of tp yields a non-constant T -periodic solution x. 



We recall once more that by the integer part E[a] of a G H, we mean the 
a e 7Z such that a < a < a + 1. For instance, if we take = 0, Corollary 2 
tells us that x exists and is non-constant provided that: 



2^ 6 ~ <1< 2^ 



or 



T G 



2tt 2tt 



(15) 
(16) 



Proof. The spectrum of A is 1^.2? + a^. The largest negative eigenvalue A is 
given by + a^, where 



— A: + aoo < < — (k„ + 1) + a c 



Hence: 

The condition 7 < —A < S now becomes: 

2tt 



aoo < 



k a - < lu - a a 



which is precisely condition (14). 



(17) 
(18) 

(19) 
□ 
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Lemma 1. Assume that H is C 2 on R 2 ™\{0} and that H"(x) is non-degenerate 
for any x ^ 0. Then any local minimizer x of ip has minimal period T. 

Proof. We know that x, or x + £ for some constant £ £ H 2 ™, is a T-periodic 
solution of the Hamiltonian system: 

x = JH'(x) . (20) 

There is no loss of generality in taking £ = 0. So ip(x) > i[)(x) for all x in 
some neighbourhood of x in IT 1 ' 2 (TR/T2Z; IR 2 ™). 

But this index is precisely the index ix{x) of the T-periodic solution x over 
the interval (0,T), as defined in Sect. 2.6. So 

ir(x) - . (21) 

Now if x has a lower period, T/k say, we would have, by Corollary 31: 

i T (x) = i kT/k (x) > ki T/k (x) + fc-l>fc-l>l. (22) 

This would contradict (21), and thus cannot happen. □ 

Notes and Comments. The results in this section are a refined version of [1] ; the 
minimality result of Proposition 14 was the first of its kind. 

To understand the nontriviality conditions, such as the one in formula (16), 
one may think of a one-parameter family xt, T £ (27rw _1 , 2irb^) of periodic 
solutions, xt(0) = Xt(T), with xt going away to infinity when T — > 2irui~ 1 , 
which is the period of the linearized system at 0. 



Table 1. This is the example table taken out of The TgXbook, p. 246 



Year 


World population 


8000 B.C. 


5,000,000 


50 A.D. 


200,000,000 


1650 A.D. 


500,000,000 


1945 A.D. 


2,300,000,000 


1980 A.D. 


4,400,000,000 



Theorem 1 (Ghoussoub-Preiss). Assume H(t,x) is {Q 1 e)-subquadratic at 
infinity for all e > 0, and T-periodic in t 

H(t, •) is convex Vi (23) 

H(-,x) is T-pcriodic Vx (24) 
H(t,x) > n(\\x\\) with n(s)s^ 1 — > oo as s — > oo (25) 
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Ve>0, 3c : H(t,x) < ^\\xf + c . (26) 

Assume also that H is C 2 , and H"(t, x) is positive definite everywhere. Then 
there is a sequence Xk, k <E IN, of kT -periodic solutions of the system 

x = JH'(t,x) (27) 

such that, for every k G IN, there is some p a € IN with: 

P>Po=>x pk ^x k . (28) 

□ 

Example 1 (External forcing ). Consider the system: 

x = JH'{x) + f(t) (29) 

where the Hamiltonian H is (0, &oo)-subquadratic, and the forcing term is a 
distribution on the circle: 

f=j t F + f with FeL 2 (M/T2Z; IR 2 ™) , (30) 
where f a := T _1 f(t)dt. For instance, 

/(*) = E s *z > ( 31 ) 

feeiN 

where S k is the Dirac mass at t = k and £ G IR 2 " is a constant, fits the pre- 
scription. This means that the system x = JH'(x) is being excited by a series 
of identical shocks at interval T. 

Definition 1. Let A^it) and i?oo(£) be symmetric operators in JR 2n , depending 
continuously on t G [0, T], such that ^(t) < B^t) for all t. 

A Borelian function H : [0, T] x JR 2n — > IR is called (A^, B^-subquadratic 



at infinity if there exists a function N(t, x) such that: 

H(t,x) = ^{A 00 {t)x,x)+N(t,x) (32) 

Vi , N(t, x) is convex with respect to x (33) 

N(t,x) > n(\\x\\) with n(s)s^ 1 -> +oo as s -> +oo (34) 

BcelR: H(t,x)<^(B 00 (t)x,x)+c Vx . (35) 



If A^t) = aool and B^t) = b^I , with < G 1R ; we shall say that 
H is (aocb^-subquadratic at infinity. As an example, the function \\x\\ a , with 
1 < a < 2, is (0, e)-subquadratic at infinity for every e > 0. Similarly, the 
Hamiltonian 

H{t,x) = ^k\\k\\ 2 + \\x\\ a (36) 
is (k, k + e)-subquadratic for every e > 0. Note that, if k < 0, it is not convex. 
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Notes and Comments. The first results on subharmonics were obtained by Ra- 
binowitz in [5], who showed the existence of infinitely many subharmonics both 
in the subquadratic and superquadratic case, with suitable growth conditions 
on H' . Again the duality approach enabled Clarke and Ekeland in [2] to treat 
the same problem in the convex-subquadratic case, with growth conditions on 
H only. 

Recently, Michalek and Tarantello (see [3] and [4]) have obtained lower bound 
on the number of subharmonics of period kT, based on symmetry considerations 
and on pinching estimates, as in Sect. 5.2 of this article. 
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1 Fixed-Period Problems: The Sublinear Case 

With this chapter, the preliminaries are over, and we begin the search for periodic 
solutions to Hamiltonian systems. All this will be done in the convex case; that 
is, we shall study the boundary-value problem 

x = JH'(t,x) 
x(0) = x(T) 

with H(t, •) a convex function of x, going to +oo when ||x|| — ► oo. 
1.1 Autonomous Systems 

In this section, we will consider the case when the Hamiltonian H(x) is au- 
tonomous. For the sake of simplicity, we shall also assume that it is C . 

We shall first consider the question of nontriviality, within the general frame- 
work of (Am, -Boo)-subquadratic Hamiltonians. In the second subsection, we shall 
look into the special case when H is (0, fe^-subquadratic, and we shall try to 
derive additional information. 

The General Case: Nontriviality. We assume that H is {A^, B^-sub- 
quadratic at infinity, for some constant symmetric matrices A^ and B^, with 
Bqc — Aoo positive definite. Set: 



7 : = smallest eigenvalue of B x — A, 



oc 



(1) 



A : = largest negative eigenvalue of J— + A, 



dt 



(2) 



Theorem 21 tells us that if A + 7 < 0, the boundary-value problem: 



x = JH'(x) 
x(0) = x(T) 



(3) 



has at least one solution x, which is found by minimizing the dual action func- 
tional: 

" ~1 



ip(u) 



-(A^u,u)+N*(-u) 



dt 



(4) 



on the range of A, which is a subspace R(A)\ with finite codimension. Here 



N(x) :=H(x)--(A oc x,x) 



is a convex function, and 



(5) 



N(x) < -((B 00 -A 00 )x,x)+c \/x 
Proposition 1. Assume H'(0) = and H(0) = 0. Set: 

S := lim inf 2N(x) \\x\\~ 2 . 

If 7 < —A < S, the solution u is non-zero: 

x{t) ^0 Vt . 



(6) 



(7) 



(8) 



Proof. Condition (7) means that, for every 5' > 5, there is some e > such that 

(9) 



\\x\\<e^N(x)< 6 -\\xf 



It is an exercise in convex analysis, into which we shall not go, to show that 
this implies that there is an i] > such that 



f\\x\\<V^N*(y)<±\\y\\ 2 



(10) 



Fig. 1. This is the caption of the figure displaying a white eagle and a white horse on 
a snow field 



Since u\ is a smooth function, we will have H/iiiiH^ < rj for h small enough, 
and inequality (10) will hold, yielding thereby: 

^i)<yi|ki|l2 + y^ll«iH 2 • (11) 
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If we choose 6' close enough to 6, the quantity ( j + jr) will be negative, and 
we end up with 

ipihm) < for h ^ small . (12) 

On the other hand, we check directly that ip(0) = 0. This shows that cannot 
be a minimizer of ■0, not even a local one. Som^O and u ^ A~ 1 (0) = 0. □ 

Corollary 1. Assume H is C 2 and {a^, boo) -subquadratic at infinity. Let £i, 
. . . , £jv be the equilibria, that is, the solutions of H' '(£) = 0. Denote by uik the 
smallest eigenvalue of H" (£k), an d se t : 



:= Min {ui, . . . ,Lu k } . 



If: 



2tt 



2tt 



T 



(13) 
(14) 



then minimization of tp yields a non-constant T -periodic solution x. 



We recall once more that by the integer part E[a] of a £ H, we mean the 
a e 2Z, such that a < a < a + 1. For instance, if we take aoo = 0, Corollary 2 
tells us that x exists and is non-constant provided that: 



or 



T , T 
— boo < 1 < ^~ 
2tt 2tt 



/ 2tt 2vr 

Te — r- 



(15) 



(16) 



Proof. The spectrum of yl is + a^. The largest negative eigenvalue A is 

given by ^fc Q + a^, where 



27T 27T., 

— + aoo < < — \ k o + 1) + Ac 



Hence: 



K = E 



T 



The condition 7 < —A < 5 now becomes: 



2tt, 



aoo < -^7^0 - aoo < w - a 



which is precisely condition (14). 



(17) 
(18) 

(19) 
□ 



Lemma 1. Assume that H is C 2 on TR 2n \{0} and that H"(x) is non-degenerate 
for any x ^ 0. Then any local minimizer x of ip has minimal period T. 
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Proof. We know that x, or x + £ for some constant £ e H n , is a T-periodic 
solution of the Hamiltonian system: 

x = JH'(x) . (20) 

There is no loss of generality in taking £ = 0. So ip(x) > ip(x) for all x in 
some neighbourhood of x in W 1 ' 2 (M/TZZ; H 2 ™). 

But this index is precisely the index ir{x) of the T-periodic solution x over 
the interval (0, T), as defined in Sect. 2.6. So 

ir(x) = . (21) 

Now if x has a lower period, T/k say, we would have, by Corollary 31: 

*t(x) = i kT /k{x) > ki T/k (x) + k - 1 > k - 1 > 1 . (22) 

This would contradict (21), and thus cannot happen. □ 

Notes and Comments. The results in this section are a refined version of 1980; 
the minimality result of Proposition 14 was the first of its kind. 

To understand the nontriviality conditions, such as the one in formula (16), 
one may think of a one-parameter family xt, T £ (27rw _1 , 27r&^ 1 ) of periodic 
solutions, xt(0) = xt(T), with xt going away to infinity when T — > 2iruj~ 1 , 
which is the period of the linearized system at 0. 



Table 1. This is the example table taken out of The TgXbook, p. 246 



Year 


World population 


8000 B.C. 


5,000,000 


50 A.D. 


200,000,000 


1650 A.D. 


500,000,000 


1945 A.D. 


2,300,000,000 


1980 A.D. 


4,400,000,000 



Theorem 1 (Ghoussoub-Preiss). Assume H(t,x) is (0, s)-subquadratic at 
infinity for all e > 0, and T-periodic in t 

H(t, •) is convex Vt (23) 

H(-,x) is T-pcriodic Vx (24) 
H(t,x) > n(\\x\\) with n(s)s^ 1 — > co as s — > co (25) 
Ve>0, 3c : H(t,x) <^\\xf + c . (26) 
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Assume also that H is C 2 , and H"(t, x) is positive definite everywhere. Then 
there is a sequence Xk, fc £ IN, of kT -periodic solutions of the system 

x = JH'(t,x) (27) 

such that, for every k £ IN, there is some p £ IN with: 

P> Po => x pk ^ x k . (28) 

□ 

Example 1 (External forcing ). Consider the system: 

x = JH'(x) + f(t) (29) 

where the Hamiltonian H is (0, &oo)-subquadratic, and the forcing term is a 
distribution on the circle: 

f=jF + f Q with F £ L 2 (M/T2Z; JR 2n ) , (30) 
where f := T _1 f(t)dt. For instance, 

/(*) = E s *z - ( 31 ) 

feeiN 

where 5k is the Dirac mass at t = k and £ £ M 2 ™ is a constant, fits the pre- 
scription. This means that the system x = JH'(x) is being excited by a series 
of identical shocks at interval T. 

Definition 1. Let A^t) and B^t) be symmetric operators in IR 2 ™, depending 
continuously on t £ [0, T], such that ^(t) < B^t) for all t. 

A Borelian function H : [0, T] x IR 2 " — > IR is called (A^, B^-subquadratic 



at infinity if there exists a function N(t, x) such that: 

H(t,x) = ^{A 00 (t)x,x)+N(t,x) (32) 

Vt , N(t, x) is convex with respect to x (33) 

N(t,x) > n(\\x\\) with n(s)s^ 1 -> +oo as s -> +oo (34) 

3c£H: H(t,x) < ^ (Boo(t)a;,a;) +c Vx . (35) 



If Aoo(t) = a^I and B^t) = b^I , with < b^ £ IR, we shall say that 
H is (aaojbooj-subquadratic at infinity. As an example, the function \\x\\ , with 
1 < a < 2, is (0, e)-subquadratic at infinity for every e > 0. Similarly, the 
Hamiltonian 

H(t,x) = h\\kf + \\x\\ a (36) 
is (k, k + e)-subquadratic for every e > 0. Note that, if k < 0, it is not convex. 
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Notes and Comments. The first results on subharmonics were obtained by Rabi- 
nowitz in 1985, who showed the existence of infinitely many subharmonics both 
in the subquadratic and superquadratic case, with suitable growth conditions 
on H 1 . Again the duality approach enabled Clarke and Ekcland in 1981 to treat 
the same problem in the convex-subquadratic case, with growth conditions on 
H only. 

Recently, Michalek and Tarantello (see Michalek, R., Tarantello, G. 1982 and 
Tarantello, G. 1983) have obtained lower bound on the number of subharmonics 
of period kT, based on symmetry considerations and on pinching estimates, as 
in Sect. 5.2 of this article. 

References 

Clarke, F., Ekeland, I.: Nonlinear oscillations and boundary-value problems for Hamil- 

tonian systems. Arch. Rat. Mech. Anal. 78, 315-333 (1982) 
Clarke, F., Ekeland, I.: Solutions pcriodiques, du periode donnee, des equations hamil- 

toniennes. Note CRAS Paris 287, 1013-1015 (1978) 
Michalek, R., Tarantello, C: Subharmonic solutions with prescribed minimal period 

for nonautonomous Hamiltonian systems. J. Diff. Eq. 72, 28-55 (1988) 
Tarantello, C: Subharmonic solutions for Hamiltonian systems via a 7Z V pseudoindex 

theory. Annali di Matematica Pura (to appear) 
Rabinowitz, P.: On subharmonic solutions of a Hamiltonian system. Comm. Pure Appl. 

Math. 33, 609-633 (1980) 



Author Index 



Abt I. 7 


Cozzika G. 9 


Ahmed T. 3 


Criegee L. 11 


Andreev V. 24 


Cvach J. 27 


Andrieu B. 27 




Arpagaus M. 34 


Dagoret S. 28 




Dainton J.B. 19 


Babaev A. 25 


Dann A.W.E. 22 


Barwolff A. 33 


Dau W.D. 16 


Ban J. 17 


Deffur E. 11 


Baranov P. 24 


Delcourt B. 26 


Barrelet E. 28 


Buono Del A. 28 


Bartel W. 11 


Devel M. 26 


Bassler U. 28 


De Roeck A. 11 


Beck H.P. 35 


Dingus P. 27 


Behrend H.-J. 11 


Dollfus C. 35 


Berger Ch. 1 


Dreis H.B. 2 


Bergstein H. 1 


Drescher A. 8 


Bernardi G. 28 


Diillmann D. 13 


Bernet R. 34 


Diinger O. 13 


Besancon M. 9 


Duhm H. 12 


Biddulph P. 22 




Binder E. 11 


Ebbinghaus R. 8 


Bischoff A. 33 


Eberle M. 12 


Blobel V. 13 


Ebert J. 32 


Borras K. 8 


Ebert T.R. 19 


Bosetti P.C. 2 


Efremenko V. 23 


Boudry V. 27 


Egli S. 35 


Brasse F. 11 


Eichenberger S. 35 


Braun U. 2 


Eichler R. 34 


Braunschweig A. 1 


Eisenhandler E. 20 


Brisson V. 26 


Ellis N.N. 3 


Biingener L. 13 


Ellison R.J. 22 


Burger J. 11 


Elsen E. 11 


Biisser F.W. 13 


Evrard E. 4 


Buniatian A. 11,37 




Buschhorn G. 25 


Favart L. 4 




Feeken D. 13 


Campbell A.J. 1 


Feist R. 11 


Carli T. 25 


Feltesse A. 9 


Charles F. 28 


Fensome I.F. 3 


Clarke D. 5 


Ferrarotto F. 31 


Clegg A.B. 18 


Flamm K. 11 


Colombo M. 8 


Flauger W. 11 


Courau A. 26 


Flieser M. 25 


Coutures Ch. 9 


Fluggc G. 2 



14 



Fomenko A. 24 
Fominykh B. 23 
Formanek J. 30 
Foster J.M. 22 
Franke G. 11 
Fretwurst E. 12 

Gabathuler E. 19 
Gamerdinger K. 25 
Garvey J. 3 
Gayler J. 11 
Gellrich A. 13 
Gennis M. 11 
Genzel H. 1 
Godfrey L. 7 
Goerlach U. 11 
Goer rich L. 6 
Gogitidze N. 24 
Goodall A.M. 19 
Gorelov I. 23 
Goritchev P. 23 
Grab C. 34 
Grassier R. 2 
Greenshaw T. 19 
Greif H. 25 
Grindhammer G. 25 

Haack J. 33 
Haidt D. 11 
Hamon O. 28 
Handschuh D. 11 
Hanlon E.M. 18 
Hapke M. 11 
Harjes J. 11 
Haydar R. 26 
Haynes W.J. 5 
Hedberg V. 21 
Heinzelmann G. 13 
Henderson R.C.W. 18 
Henschel H. 33 
Herynek I. 29 
Hildesheim W. 11 
Hill P. 11 
Hilton CD. 22 
Hoeger K.C. 22 
Huet Ph. 4 
Hufnagel H. 14 
Huot N. 28 

Itterbeck H. 1 



Jabiol M.-A. 9 
Jacholkowska A. 26 
Jacobsson C. 21 
Jansen T. 11 
Jonsson L. 21 
Johannscn A. 13 
Johnson D.P. 4 
Jung H. 2 

Kalmus P.I. P. 20 
Kasarian S. 11 
Kaschowitz R. 2 
Kathage U. 16 
Kaufmann H. 33 
Kenyon I.R. 3 
Kermiche S. 26 
Kiesling C. 25 
Klein M. 33 
Kleinwort C. 13 
Knies G. 11 
Ko W. 7 
Kohler T. 1 
Kolanoski H. 8 
Kole F. 7 
Kolya S.D. 22 
Korbel V. 11 
Korn M. 8 
Kostka P. 33 
Kotclnikov S.K. 24 
Krehbiel H. 11 
Kriicker D. 2 
Kriiger U. 11 
Kubenka J. P. 25 
Kuhlen M. 25 
Kurca T. 17 
Kurzhofer J. 8 
Kuznik B. 32 

Lamarchc F. 27 
Lander R. 7 
Landon M.P.J. 20 
Lange W. 33 
Lanius P. 25 
Laporte J.F. 9 
Lebedev A. 24 
Leuschner A. 11 
Levonian S. 11,24 
Lewin D. 11 
Ley Ch. 2 
Lindner A. 8 



Lindstrom G. 12 


Patel G.D. 19 


Linsel F. 11 


Peppel E. 11 


Lipinski J. 13 


Phillips H.T. 3 


Loch P. 11 


Phillips J. P. 22 


Lohmander H. 21 


Pichler Ch. 12 


Lopez G.C. 20 


Pilgram W. 2 




Pitzl D. 34 


Magnussen N. 32 


Prell S. 11 


Mani S. 7 


Prosi R. 11 


Marage P. 4 




Marshall R. 22 


Radel G. 11 


Martens J. 32 


Raupach F. 1 


Martin A.® 19 


Rauschnabel K. 8 


Martyn H.-U. 1 


Reinshagen S. 11 


Martyniak J. 6 


Ribarics P. 25 


Masson S. 2 


Riech V. 12 


Mavroidis A. 20 


Riedlberger J. 34 


McMahon S.J. 19 


Rietz M. 2 


Mehta A. 22 


Robertson S.M. 3 


Meier K. 15 


Robmann P. 35 


Mercer D. 22 


Roosen R. 4 


Merz T. 11 


Royon C. 9 


Meyer C.A. 35 


Rudowicz M. 25 


Meyer H. 32 


Rusakov S. 24 


Meyer J. 11 


Rybicki K. 6 


Mikocki S. 6,26 




Milone V. 31 


Sahlmann N. 2 


Moreau F. 27 


Sanchez E. 25 


Moreels J. 4 


Savitsky M. 11 


Morris J.V. 5 


Schacht P. 25 


Miiller K. 35 


Schleper P. 14 


Murray S.A. 22 


von Schlippe W. 5 




Schmidt D. 32 


Nagovizin V. 23 


Schmitz W. 2 


Naroska B. 13 


Schoning A. 11 


Naumann Th. 33 


Schroder V. 11 


Newton D. 18 


Schulz M. 11 


Neyret D. 28 


Schwab B. 14 


Nguyen A. 28 


Schwind A. 33 


Niebergall F. 13 


Seehausen U. 13 


Nisius R. 1 


Sell R. 11 


Nowak G. 6 


Semenov A. 23 


Nyberg M. 21 


Shekelyan V. 23 




Shooshtari H. 25 


Oberlack H. 25 


Shtarkov L.N. 24 


Obrock U. 8 


Siegmon G. 16 


Olsson J.E. 11 


Siewert U. 16 


Ould-Saada F. 13 


Skillicorn I.O. 10 




Smirnov P. 24 


Pascaud C. 26 


Smith J.R. 7 



16 



Smolik L. 11 
Spitzer H. 13 
Staroba P. 29 
Steenbock M. 13 
Steffen P. 11 
Stella B. 31 
Stephens K. 22 
Stosslein U. 33 
Strachota J. 11 
Straumann U. 35 
Struczinski W. 2 

Taylor R.E. 36,26 
Tchernyshov V. 23 
Thiebaux C. 27 
Thompson G. 20 
Truol P. 35 
Turnau J. 6 

Urban L. 25 
Usik A. 24 

Valkarova A. 30 
Vallee C. 28 
Van Esch P. 4 
Vartapetian A. 11 



Vazdik Y. 24 
Verrecchia P. 9 
Vick R. 13 
Vogel E. 1 

Wacker K. 8 
Walther A. 8 
Weber G. 13 
Wegner A. 11 
Wellisch H.P 25 
West L.R. 3 
Willard S. 7 
Winde M. 33 
Winter G.-G. 11 
Wolff Th. 34 
Wright A.E. 22 
Wulff N. 11 

Yiou T.P. 28 

Zacek J. 30 
Zeitnitz C. 12 
Ziaeepour H. 26 
Zimmer M. 11 
Zimmermann W. 



Subject Index 



Absorption 327 

Absorption of radiation 289-292, 299, 
300 

Actinides 244 

Aharonov-Bohm effect 142-146 
Angular momentum 101-112 

- algebraic treatment 391-396 
Angular momentum addition 185-193 
Angular momentum commutation rela- 
tions 101 

Angular momentum quantization 9-10, 
104-106 

Angular momentum states 107, 321, 
391-396 
Antiquark 83 
Q-rays 101-103 

Atomic theory 8-10, 219-249, 327 
Average value 

(see also Expectation value) 15-16, 25, 34, 
37, 357 

Baker-Hausdorff formula 23 
Balmcr formula 8 
Balmer series 125 
Baryon 220, 224 
Basis 98 

Basis system 164, 376 

Bell inequality 379-381, 382 

Bessel functions 201, 313, 337 

- spherical 304-306, 309, 313-314, 322 
Bound state 73-74, 78-79, 116-118, 202, 
267, 273, 306, 348, 351 

Boundary conditions 59, 70 
Bra 159 

Breit-Wigner formula 80, 84, 332 



Brillouin-Wigner perturbation theory 
203 

Cathode rays 8 
Causality 357-359 
Center-of-mass frame 232, 274, 338 
Central potential 113-135, 303-314 
Centrifugal potential 115-116, 323 
Characteristic function 33 
Clebsch-Gordan coefficients 191-193 
Cold emission 88 
Combination principle, Ritz's 124 
Commutation relations 27, 44, 353, 391 
Commutator 21-22, 27, 44, 344 
Compatibility of measurements 99 
Complete orthonormal set 31,40,160, 
360 

Complete orthonormal system, see 
Complete orthonormal set 
Complete set of observables, see Complete 
set of operators 

Eigenfunction 34, 46, 344-346 
- radial 321 

- calculation 322-324 
EPR argument 377-378 
Exchange term 228,231,237,241,268, 
272 

/-sum rule 302 
Fermi energy 223 

H^" molecule 26 
Half-life 65 
Holzwarth energies 68 



